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TECHNICAL PAPER 


SPACECRAFT FLIGHT CONTROL SYSTEM DESIGN SELECTION PROCESS 
FOR A GEOSTATIONARY COMMUNICATION SATELLITE 

I. INTRODUCTION 


The Earth’s first artificial satellite, Sputnik I, slowly tumbled in orbit. The first U.S. satellite, 
Explorer I, also tumbled out of control. 1 In 1962 Alouette /, and in 1964 Explorer XX, both suffered 
rapid spin decay due to solar torque. In 1969, the Applications Technology Satellite 5 (ATS5) quickly 
began to tumble and went unstable due to energy dissipation in the heat pipe, and entered a flat 
spin. 2 Today, as we launch the Mars observer and the Cassini spacecraft, and with missions 
demanding increasing life expectancies, satellite stability and control has become a much higher 
priority. 

In order to select an attitude control system for a spacecraft (S/C) and to investigate the 
stability requirements of that S/C, first the anticipated disturbance torques must be identified and 
quantified. Next, appropriate attitude control system options to counteract these torques must be 
reviewed according to the specifications of the S/C mission, and the dynamic stability of the S/C 
verified. To see the process, a geostationary communication satellite with a life expectancy of 10 to 
14 years is used as an example. 


II. ASSESSMENT OF DISTURBANCE TORQUES 


The anticipated disturbance torques consist of environmental torques and torques internal to 
the S/C. 

Aerodynamic . The interaction of the upper atmosphere with a S/C’s surface produces an 
aerodynamic torque about its center of mass. For S/C below 400 km, 3 the aerodynamic torque is the 
dominant environmental disturbance torque. It is customary to neglect aerodynamic effects above 
500 km. In geostationary orbit (GEO) where the altitude is 35,786 km, this torque can be safely 
neglected. 

Ma gnetic . Magnetic disturbance torques result from the interaction between the S/C’s 
residual magnetic field and the geomagnetic field. The primary sources of magnetic disturbance 
torques are: S/C magnetic moments, eddy currents, and hysteresis. Of these, the S/C’s magnetic 
moment is usually the dominant source of disturbance torque. This disturbance torque is a function of 
the geocentric magnetic flux density, and is of significant consideration below 1,500 km. Since the 
flux density decreases with altitude, the magnetic torque can be neglected for an S/C in GEO. 

Gravity Gradient . Any nonsymmetrical S/C in orbit is subject to a gravitational torque 
because of the variation in the Earth’s gravitational force over the S/C body. A spherical or non- 
spherical Earth model can be used. The noncircularity of the equatorial cross section of the Earth, as 
expressed by the tesseral terms of the Earth’s gravitational potential, causes a longitude drift rate 
which is a function of longitude over which the satellite is stationed. For a satellite whose longitude 
is maintained within a narrow band by station-keeping maneuvers, this tesseral acceleration can be 



assumed to be constant. Gravity gradient torque is normal to the local vertical and vanishes for a 
spherically symmetric S/C. The gravity-gradient disturbance torque for a low Earth orbit (LEO), 
weak at best, is over 200 times weaker at GEO, and can be neglected. 

Collision . Longitudinal separations of neighboring S/C in GEO are regulated by the 
International Telecommunication Union to restrict both intersatellite and intersystem interference. 
Studies are being conducted to investigate the possibilities of sharing an allocated longitudinal slot 
by more than one satellite. 4 Current longitudinal band separation is ±0.1° (74 km). This sharing of 
assigned longitudinal slots will increase chances of collision and demand more accurate S/C control 
systems. A communication satellite in a shared slot would need increased accuracy. More stringent 
specifications in terms of jitter and drift rate and higher stability are covered in the literature. 5 
Collision torque is neglected in this assessment. 

Debris. Micrometeorite . Encounters with large particles can degrade S/C surface properties in 
time, can deteriorate thermal coatings and solar cells, and small high speed particles can cause 
cratering. However, these encounters usually do not produce significant disturbance torques. For 
this assessment, the disturbance torque from debris and micrometeorites is neglected. 

Solar Radiation Torque . Because the solar radiation varies as the inverse square of the 
distance from the Sun, the solar radiation pressure is essentially altitude independent for S/C in 
Earth orbit. The major factors determining the radiation torque on a S/C are the intensity and spectral 
distribution of the incident radiation, the geometry of the surface and its optical properties, and the 
orientation of the Sun vector relative to the S/C. Torque due to solar radiation pressure, T S p, is highly 
dependent on the type of surface being illuminated, and also on the S/C geometry. An ideal surface is 
either transparent, absorbent, or a reflector, but most S/C surfaces are a combination of all three. 
Reflectors are diffuse or specular. In general, solar arrays are absorbers and the S/C body is a 
reflector. 

The torque due to solar radiation can be estimated by 

T sp = Ps As Ls(l+q) cos i , 


where 


p s =4.617xl0 -6 N/m 2 = solar constant 
As = area of S/C surface 

Ls = center-of-pressure to center-of-mass offset 
i = angle of incidence of Sun 
q = 0.6 = reflectance factor . 6 

Internal Torques . Internal torques are defined as torques exerted on the main body of a S/C 
by internal moving parts such as its reaction wheels, flexible booms, or partially filled propellant 
tanks. In the absence of external torques, the total angular momentum remains constant. However, 
internal torques can alter the system’s kinetic energy and redistribute the S/C’s angular momentum 
among its component parts in ways which can change its dynamic characteristics. 3 For example, in a 
spinning S/C, angular momentum can be transferred from the nominal spin axis to another principal 
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axis resulting in nutation, uncontrolled tumbling, or flat spin. These undesirable results are often 
best’ countered by attitude-stabilization systems based on other internally generated torques sue 
a S ^Lacflon control systems (RCS’s), nutation dampers, reaction wheels, and other movable-mass 

stabilizing mechanisms. 


The most important anticipated S/C disturbance torques have been assessed. At the design 
altitude the quantified dominant disturbance torque has been determined to be the solar iadia 
torque. ’ Figure 1 shows the relative magnitudes of the considered environmental disturbance 

torques. 1 



Figure 1. Environmental disturbance torque comparisons. 


III. CONTROL TORQUE OPTIONS 


Snin Stabilization . This is the most common passive control technique in which the entire S/C 
is rotated so that the momentum vector remains approximately fixed in inertial spa«. 
method requires the use of an active control system to periodically adjust the S/C alt tude and spin 
rale ancTto counteract disturbance torques. RCS’s are commonly used A — of wo 
reorientation thrusters and two spin rate control thrusters are required^ n " 

also requires nutation damping where the nutation is caused by an unbalanced " ^ ^ ^ f 

the S/C structure. The S/C normally spins about its major principal axis for stability and is caiiea 
singk st^n conflol. Spin stabilization is a simple and effective control technique and requires no 
moving parts, howevei? it is limited to S/C for which the spin itself does not inhibit the S/C function. 
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Dual Spin Control . The S/C has a rotating wheel and a platform, or consists of two rotating 
components spurning at different rates. This is a passive control system, but it reSes the o J an 

active control system to adjust the S/C attitude and spin rate and to counteract disturbance torques 

n 1S , erf 6 Contro system a,so requires nutation damping due to imbalances in the S/C and elasticity. 
- pin control operates on the same principle as a single-spin, but provides platforms for both 

p0 . inting instraments - Wd h a two-component S/C, additional complexity arises because 
of the need for bearings and support structures separating the two components. 

gravity Gradient, The differential gravitational forces acting on an asymmetric S/C force the 
m,no r ax,s to be perpendicular to the gravitational equipotemial. The pointing accuracy ia LiSfy r 

nd helT onT nro ‘ S/C Where °" e is ^ 

and best on a prolate body, causing the minor axis to align along the nadir vector. To obtain the 

“ in — S ° inerUa ’ b °° mS are 0ften de P loyed along the minor axis. The gravity- 
g ent control torque causes the S/C to oscillate or librate about the pitch axis, and a passive 
amper is generally required to minimize the amplitude of the oscillation. Gravity-gradient control 

STachTpr 6 n ° ?h° VinS PartS °! her than in S ° me Cases extendab,e booms which can be 230-m 
ong each_ Because the gravity-gradient torque decreases as the inverse cube of the vector distance 
from the Earth, these control systems are best for LEO. mstance 

K utation Damping A nutation damper aligns the spin axis with the angular momentum vector 
by dissipating the excess kinetic energy and moves the S/C by removing the nutational motion A 
^mentum control system will tend to exhibit nutation when within the RCS deadbands The 
amplitude of this nutation is dependent on the minimum impulse bit for the RCS. 7 A new method of 

D^vi^d S "° l reqUire m ° memUm deSaluration ha s been proposed by 

, . Pa&> lv e Versus Active Control . The most common active S/C control systems involve mass 
expulsion devices, such as gas jets or ion thrusters; momentum wheels; and electromagnetic coils 

lf!-h, nCe 7 H PreSent , S 3 Vanety ° f momentum control systems. Greater chicle con^S^SSSlS^ 
possible, and control torques are larger with active control. y 

Magnetic Coil . These control systems use magnetic torquing for control. Electromagnets and 
iron-core magnets are used for both stabilization and maneuvering Electromagnet!^ 
vary the control coils’ polarity and direction to match the Earth’s magnetic fiefd to produce a torque 
to change attitude. In most systems, a set of three mutually perpendicular coils are used 3 Some S/C 

coTcX/nTl'T ° f n,ag " eliC COi ' conm>l and RCS 10 a momentum whed system Znefc 

control systems require no moving parts, complex hardware, or expendables, however they are 
slow in maneuvering and also, only appropriate for S/C at altitudes less than GEO. ’ V 

ct .... Momentum Control . Momentum wheel control systems can be used for S/C maneuvering and 
stabilization They maintain attitude by momentum exchange between the S/C and the wheel S As a 
disturbance torque acts on the S/C along one axis, the wheel reacts, absorbing the torque and 

Mom^ntlinf wtf f md t e - The wheel spin rate increases or decreases to maintain alonsmtZltude 
omentum wheel control systems are particularly suited for attitude control in the presence of cvclic 
torques or random torques. Secular torques cause the wheel speed to either increase Sr Secrea S 
monotonically until the wheel speed is outside the operational constraints. An active control system 
must then be used to restore the momentum wheel speed to its nominal operating value MomSm 
wheel control systems can have wheels on one, two, or three axes. Figure 2 shows several possible 
configurMions. A two-wheel system for an Earth oriented S/C normally has one wheel along the 
pitch axis for pitch control and the other on either the roll or the yaw axis for roll/yaw control 
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Figure 2. Various momentum control systems. 
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Gimbaled momentum wheels change both momentum magnitude and direction. The double-gimbaled 
momentum wheel permits a continuous rotation of a S/C if the nominal direction of the momentum 
vector and the S/C’s rotation axis are aligned. The three-axis control provided by the gimbaled 
momentum wheel control system is well suited for S/C that must be aligned periodically to different 
ground stations. Momentum control systems have proven performance and relative simplicity, permit 
high accuracies, operate more smoothly than mass expulsion systems, and produce no exhaust 
contaminants. They can be mixed and matched to suit specific requirements, but do need auxiliary 
torque to desaturate the momentum wheels. 

Bias Momentum . In this type of control system, one or more wheels provide a bias or nonzero 
angular momentum. Gyroscopic stiffness is provided against periodic disturbance torques, and 
variable momentum storage capability is provided about its spin axis. Small, high-speed wheels are 
usually preferred due to their low size and weight. Advantages of this control system are short-term 
stability against disturbance torques similar to spin stabilization, roll/yaw coupling that permits yaw 
angle stabilization without a yaw sensor, and a momentum wheel that may be used as an actuator 
for pitch angle control. Thus, a momentum bias control system can provide three-axis control with 
less instrumentation than a three-axis reaction wheel. 

Control Moment Gvro (CMG) . This control system can be a single or double gimbaled wheel 
spinning at a constant rate. Single-gimbal CMG’s are used for high torque maneuvering. Double- 
gimbal CMG’s are used to absorb cyclic disturbance torques. This method of S/C control is used on 
the MIR Space Station and Space Station Freedom, and is most suitable where large control torques 
are required. 

Reaction Wheels. Reaction torques are generated about the spin axis by accelerating or 
decelerating a wheel from its nominal spin rate of zero. The reaction torque causes the S/C to rotate 
in the opposite direction. Three orthogonal reaction wheels can provide three-axis control and high 
pointing accuracy. Reaction wheel control systems are well suited for cyclic torque absorption and for 
low torque momentum transfer during reorientation maneuvers . 7 Because the disturbance torques in 
GEO are small, it is possible to use small reaction wheels to absorb them with an active control 
system to maintain three-axis stability. Advantages of three-axis reaction wheel control systems 
are capability of continuous high-accuracy pointing control, large angle slewing maneuvers without 
fuel consumption, and compensation for cyclic torques without fuel consumption. Configurations of 
more than three reaction wheels can be used providing backups. In this case, a steering law is 
required to distribute the momentum between the wheels when all are operating simultaneously. For 
example, the Hubble space telescope uses a four-wheel reaction control system. 

Three-Axis Control . Attitude acquisition can take up to a few weeks and involves the 
maneuvering necessary to reorient and reconfigure the S/C from tip-off (separation from the launch 
vehicle) to its mission operation mode. The control problems unique to attitude acquisition include 
the deployment of extendable booms, antennas, solar panels, and flight control surfaces. A new 
method of attitude acquisition for a communication satellite has been developed using quaternion 
feedback to autonomously despin the S/C . 9 For a short mission, three-axis control can be achieved 
entirely by an RCS. Thrusters are usually fired in pairs to minimize translational motion. RCS 
thrusters are efficient in maneuvering the S/C, are simple to operate, and are not limited to a specific 
environment. However, they are more costly, require more hardware, and are limited by fuel onboard. 
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For a mission life of 10 to 14 years, a control system composed entirely of RCS is not 
recommended. Various combinations of momentum control systems can be used to provide three- 
axis control. A three-axis momentum control system has momentum wheels along all three axes and 
may have six or more wheels along nonorthogonal axes. Figure 3 lists various three-axis control 
system configurations. Reference 14 incorporated a hybrid of these. 


SINGLE MOMENTUM WHEEL PITCH MOMENTUM WHEEL/YAW REACTION WHEEL 

PITCH MOMENTUM WHEEL/THRUSTER DOUBLE-GIMBAL MOMENTUM WHEEL 

SINGLE-GIMBAL MOMENTUM WHEEL THREE REACTION WHEELS 

CANTED SCAN WHEEL MOMENTUM BIAS 


Figure 3. Three-axis control configurations. 


IV. FLIGHT CONTROL SYSTEM DESIGN 


After assessing the disturbance torque environment, and the S/C mission, and reviewing the 
control systems available, the selected flight control system design is: 

An active, three-axix control, double-gimbaled, bias momentum control system with 

a. RCS 

or 

b. Solar flight control surfaces 
for desaturation. 

The three-axis control provided by the gimbaled momentum wheel control system is well 
suited for satellites that must be aligned periodically to different ground stations and can be powered 
electrically from solar cells or onboard generators. Options A or B can be used for momentum 
desaturation of the wheel. In sizing the RCS, the duration of the desaturation impulses is a function 
of the amount of excess momentum. This is typically about 1 percent of the nominal wheel momen- 
tum. Efforts to reduce thruster usage on GEO satellites in order to reduce fuel consumption are well 
documented. 10 

Option B uses flight control surfaces which put the dominant disturbance torque, the solar 
disturbance torque, to advantage to provide vehicle control and wheel desaturation, instead of using 
an RCS. Figure 4 shows a recommended flight control system, and figure 5 shows a typical 
configuration of a S/C using flight control surfaces. This particular configuration has eight flight 
control surfaces for attitude stabilization and control. 11 
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Figure 4. Selected flight control system. 

SUN SENSOR 



Figure 5. Option B using flight control surfaces for desaturation 




Euler’s moment equations for this S/C control system design are, 

T=h = h b +co xh ; h = h v +h w ; h v = I x co x i+I y (o y j+l z O) z k 
T = [I x d) x -S (sin 8 sin y)h w + y (cos 8 cos y)h w +( cos 8 sin y)h w 

+ co ( col-K sin 8)-(0 z {(0 l y -h w cos 8 cos y)] i 

y z y 

+[I y 0) y +8 (sin 8 cos ?)h w +Y (cos 8 sin y)h w -(cos < 5 cos y>h w 
+co (co I x +h w cos 8 sin y)-(o x {coI z -h w sin 8)]j 
+[l,6) z -8 (cos 5 )/i M ,-(sin 8)h w +0) x (a) y I y -h w cos 5 cos y) 
-O) ( 0 ) I x +h w cos 5 sin y)]£ , (ref. 2 ) 

where 

7 = the total disturbance torque on the S/C 
h = the total angular momentum 

h b = the rate of change of angular momentum with respect to the body 
ta = angular velocity 


h v = angular momentum of vehicle 


h w = angular momentum of wheel 

8 = roll gimbal angle 

y- yaw gimbal angle. 

Wheel momentum components are, 

h wx = (cos 8 sin y)h w ; h ny = -(cos < 5 cos y)h w ; h wz = -(sin 8 )h w . 

These equations are highly nonlinear, and Lyapunov’s direct method is still one of the most 
powerful techniques available today for stability analysis of nonlinear systems. 

For small gimbal deflections, small deviations from nominal wheel momentum, and 


h n »max[I x co o ,I y a) 0 ,I z co f ) , 

and l x = l z , the above becomes three equations in pitch 6 , roll <p, and yaw yr. Setting two of the three 
principal moments of inertia equal avoids the asymmetric case, where the solution of Euler s moment 
equations cannot be written in terms of trigonometric functions. In the asymmetric case the solu 
involves the Jacobian elliptic functions (see references 15 , 16 , and 7 ). 
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Pitch: T y = I y 6+h y , 


Roll: T x - I x <p +Q) o h n <p+h n \f/+h x ~ah z , 
Yaw: T z = I z \ l /+ Q) h n¥ -h n <i) +h z -coJi x , 


where h y ,h x ,h z are control functions 
h„ = nominal wheel momentum 
0 ) o = GEO orbital rate = 7.28xl0- 5 rad/s 


Solar pressure torques (t = 0 at 6 a.m./p.m.) 


T x = 2x1 0“ 5 (1-2 sin co 0 t) Nm 
T y = 1(H (cos ( 0 o t ) Nm 
T z = -5x1 0~ 5 (cos (D 0 t ) Nm 
lx = I z = 3,000 Nms 2 
I y = 660 Nms 2 
S/C mass = 1 ,000 kg 
Pitch and roll accuracy = 0.05° 

Yaw accuracy = 0.40°. 

PHgh Control . To provide the required pitch damping, a pitch control law is selected 


h y = Kp(T p 0+d ) , 

where K P = pitch autopilot gain in Nm/rad and t= time constant in seconds. Therefore, 

T y = I y e+K P T p d+K P e => 

T + 





■V <t,s + n fe 


PSEUDORATE MODULATOR 

Since the sensors employed can provide only direct measurement of angles and damping requires 
taMstep a " S raleS ' a pseudorate modulator using a Schmitt trigger with hysteresis is^sed 
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1 


0(S) 

T(S) I y S 2 +K P r p S+Kp ^ 

(Rate gyros are not appropriate here because they cannot detect very low rate while satisfying the 
long-life requirement.) 



where (0 P = pitch motion natural frequency and fj> = pitch damping. 

Selection of pitch autopilot gain Kp is based on steady-state error and response time. 
Steady-state errors must be well below the pitch accuracy limit. The maximum magnitude of solar 
pressure torque is expected to be 1(H Nm about the pitch axis. Assuming x P « orbital period, and 
the critical damping Cp = 1. the maximum steady-state error is estimated through the final value 
theorem (S = 0) as 


-4 

d ss = - T rad => 10 ~ < 0.000873 


T 

_Z 


Kp 

^pit ch > 0-1145 Nm/rad 


Since, 


Tp=1 \/k. < *pki. =01145Nm/rad=: ' T P M =152seconds - 


Fast response times are desirable, and larger values of Kp will result in faster response time. 
Choose 



which will give a response time of 80 seconds. 

Roll/Yaw Control . Control torques are produced through gimbal deflections. Damping can be 
provided through selection of an actuator control law. This will also yield a steady-state yaw error 
whose magnitude is a function of h„. 

Active roll/yaw control is investigated using the above roll and yaw equations. Roll response 
to disturbance torques should be fast and well damped. This will also minimize coupling of roll errors 
into yaw. The roll control law used is, 

M xc = h x -o)h z = Kt<(> +K(f>- 0 )h n <t> , 


where M xc is the commanded roll control torque. 
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The ideal yaw response should be well damped and decoupled from roll. The yaw control law 
used is, 

M zc = h z +co o h x -hj = -kK(r<p +<p) 

where M zc = commanded yaw control torque and k = yaw-to-roll gain ratio. 

When orbit rate decoupled momentum commands are used, the above control laws become: 


M xc = h X(r o)h zd = K(t<P +<p) , 

M zc = h zd +coh xd = -kK(t<p +<p) . 

Using a Laplace transformation on this set of equations and solving the developed transfer function, 



The high frequency roots (ft) r oii> Croii) characteristic of roll dynamics, while the low frequency 
roots (6>yaw» Cyaw ) dominate the yaw motion during yaw error correction. A large value of h n will 
result in fast yaw corrections with respect to orbit period. 


The roll autopilot gain K R is selected to limit steady-state roll error caused by a constant roll 

torque and to provide a fast response. The steady-state roll error produced by a constant roll torque 
is 


Based on the maximum roll torque resulting from solar pressure and the allowable roll error of 0 05° 
this yields 

A'ron >0.07 nm/rad . 


^roii - 0.07 Nm/rad =» r ro n = 414 seconds , 


where 

Choose, 


which will give a response time in roll of 80 seconds. 

t = 80 s => o = 0.025 rad/s . 



^roii = 1-875 Nm/rad 
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The solutions in yaw 


control from above are: 


For 


The above satisfies 



h n = 200 Nms , k = 0.0033 = 



*ro!l 


Kyw = 0.06188 Nm/rad . 


^roll T roll^ > 

1.875 (80) (3,000) > (0.033) (200) (3,000) 
150 > 6.6 => 


h = 200 Nms . 


Using M = 3.2 h A the mass of the wheel housing and electronics is 27 kg. Also : 

T z +kT x 5x10~ 5 +(0-Q33) 6x10 5 
C0 o K ~ 7.28xl0 -5 (200) 


yi/ ss = 0.0036 rad = 0.20° , 


which is well below the yaw accuracy requirement of 0.40°. 


although only three thrusters would be needed if thrust vectoring were emp oy 
momentum saturation could be relieved with the flight control surfaces. 
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V. S/C DYNAMIC STABILITY 


In a S/C with a passive control system, such as a simple spin system an unstable design will 
precess, then wobble, and finally tumble out of control due to die 

wobble can generate irregular signals from a satellite. It has been well documented that dissinative 

ar^stahHiV^^n^ 103 ^^ P rolate , S/C make them tumble out of control, while dynamically oblate S/C 
are stabilized spinning about their axis of symmetry. If the nominal S/C SDin axis is mhpr than th* 

major principal axis, energy dissipation will result in an increase in nutation which will finally 
becorae rotatton about an axis perpendicular to the original axis and known as a flaTspin Recover 
*P ln . re r res ^orientation of 90” and has been approached in the Hterature ustaf 

S N /cts na ' S ' abi ' ity ^ be “ S ‘“ died “ si "S ^“seco“d 8 

s ,ng 

more parts thm^n Ss^ttSTd,” m “ de P ° SSib ' e by consl ™ ctin S the S/C from two or 

A communication satellite attitude acquisition system can have reduced nrnneiiam 
=n na ^T^ WOrS by SCalin ^ the feedback ™ntrol system gaim; V proporticml to°the S/C 
COUROHS 1 aS* IZ :T'l’ Lyapun0v Stability ana,ysis ca " «» to sho P r,w"tpe of g/fn 
been correctly selected, momentum exchange ^ IcelJo,, of 

nom° dl ^ momentum adjustment is necessary due to secular disturbance torques for on-slation 
°™ ina • operation. During the attitude acquisition sequence, the S/C can be spin-stabilized in 
apogee transfer, before the momentum wheel control system is activatpd Thp Hpc-c . 

° f Sains and wheel momentum 


VI. CONCLUSIONS 


torques' has ieen'conducS “d ItedSS 

(2) Appropriate S/C flight control systems for a communication satellite have been reviewed 

xtsssjsi™ sss sr ■*" * * “ 
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